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ABSTRACT 

We develop a simple mean field approach to the transport of a passive scalar 
for which the fundamental equation is a second order differential equation 



^ I in the transported quantity, not a first order equation. Triple correlations 

CT) ' are included, as they must be for any realistic description of turbulence. No 



correlation time enters the theory, only an eddy turnover time. The approach 



CN I is simpler than standard approaches which incorporate triple correlations, but 

^ ' more realistic than Gaussian or short correlation time closures which do not. A 

O ■ similar approach has proven useful in magnetohydrodynamics. 

^ ■ 

PACS codes: 47.27.Qb; 92.10.Lq; 94.10.Lf; 95.30.Qd;96.20.Br;83.50.Xa 

o : 

^ , Transport of a passive scalar in a turbulent flow is a classic problem in the physics of 

^ ! fluids [e.g. 1-10]. Most recent formal work has focused on spectral scahngs, non-Gaussian 

statistics, and high order correlations of a passive scalar fluctuations [e.g. 9, 10] rather 
^ i than on the dynamical theory of the evolution of a mean scalar fleld [1,11-14]. However, 
^ [ understanding how a mean scalar such as average concentration evolves in a turbulent flow 
is important for a broad range of problems such as the transport of pollutants and cloud 
droplets in atmospheric turbulence [13], dust transport in planet forming disks [14], and 
particle transport in industrial flows [15]. 

A proper mean field theory requires a closure which allows the evolution of the mean 
scalar to be represented by a finite set of solvable realistic equations. If a closure can be 
tested and proven useful for a the problem of mean scalar transport, further application and 
testing in fluid dynamics and magnetohydrodynamics (MHD) is motivated. For problems 
amenable to a mean fleld treatment, a simpler closure than that needed for a full spectral 
treatment may be sufficient. Here we derive a simple practical approach to the Eulerian 
transport of a mean scalar in a turbulent flow that differs from standard treatments [e.g. 
1-8]. but whose analogy has proven effective in mean-fleld magnetic dynamo theory [16]. We 
also discuss the physical motivation for the approach and the relation to previous closures. 
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To study mean scalar transport, we start with the basic conservation equation [1, 6, 7] 
for the concentration of a passive contaminant, C. For incompressible flow, we have 

a^c- -V vc-v- vc, (1) 

where C = c + C, and c and C are the fluctuating and mean components of C respectively. 
The velocity has only a fluctuating component v whose mean vanishes. Averaging (1) gives 

dtC = -v^. (2) 

To proceed, we need to derive an equation that expresses the right side of (2) in terms of 
C. This requires the evolution equation for c. Subtracting (2) from (1) we obtain 

a^c = -V • VC - V • Vc + aT^Vc. (3) 

We are now faced with the need for a closure: If we re-write the term on the right of (2) as 
/ V ■ \/dtc{t')dt' , then using (3), a triple correlation of the form / v(t) ■ V(v(t') ■ Vc{t'))dt' 
arises. In this approach, triple correlations are sometimes ignored by one of the following 
arguments: (i) c « C, (ii) Gaussian statistics [e.g. see 5, 7] for which correlations of odd 
number of fluctuating quantities vanish, or (iii) the (5-function-correlated in time [3, 4, 17] 
(or flnite but short correlation time [11, 18]) approximation. These short correlation time 
approximations relate to the fact that upon one more iteration of the triple correlation 
by substituting c(t') = / dt"c{t")dt" , and using (3) again, the triple correlations can be 
shown to contribute a term of order Tc/rgd times that of the double correlation that 
arises from plugging the second term of (3) back into (2), where is the correlation time 
and Ted = 1/^2/^2 is an eddy turnover time associated with speed V2 = /^/S, and 
wavenumber, /c2, of the dominant turbulent scale. If Tc/Ted < 1 then triple correlations for 
quantities at two different times might be argued to be small. Using one of the above three 
approximations, and after assuming isotropy, (2) would become a diffusion equation 

dtC = (1/3) j w{t) ■ w{t')dt'V^C ~ (Te/3)^V2C = pV^C, (4) 

where the similarity follows from replacing the time integral with multiplication by Tc and 
13 is the diffusion coefficient. The solution for C is then 

C = CoE^v[-{t/red){T,/Ted){kl/kl)i (5) 

where ki is the wavenumber associated with C. Eq. (5) reduces to CoExp[— (t/3red)(/ci//i;|)] 
as Tc — > Ted. Because we focus on Eulerian transport, obtaining (4) required the signiflcant 
approximations discussed above. (For Lagrangian transport, (4) would follow more simply 
[e.g. 1-2,6-8].) But each of the arguments that led to neglecting Eulerian triple correlations 



-3- 



is poorly justified. First, c « C cannot be guaranteed at all times. Second, although 
Gaussian statistics of fluctuating quantities enforce triple correlations to vanish exactly 
for any correlation time, finite triple correlations arc required in order for there to be any 
nonlinear turbulent cascade [see 5, 6]. Third, Tc/red < 1 is not guaranteed, thus invalidating 
a universally justified expansion in this quantity. 

To develop an improved Eulcrian approach which avoids the above weaknesses, we 
re-start with (2) and instead take the time-derivative of dtC to obtain 

Si^C = -dtV^c = -9tv • Vc - V • Vdtc. (6) 

Using (3) in the last term of (6) we have 

—V • Vdtc = ujvjdidjC + VidiVjdjC + v • V(v • V)c, (7) 

where the last term in (3) does not contribute when placed in the average. The second 
term on the right vanishes for homogeneous incompressible turbulence. The first term on 
the right, after isotropizing, gives the diffusion term, and the last term is a triple correlation 
which we denote by Tc. Plugging the result back into (6) we then have 

d^C = -dtV^ = -dtv ■ Vc + Tc + ^v^V^a (8) 

We must now use the equation of motion 

d^v = - V Vv - Vp, (9) 

(where p is the pressure, and we assume incompressibility and set density p — 1) for the 
first term on the right of (8). Noting that V^p = —{diVj){djVi) for incompressible flows, we 

see that p and thus Vp depend on quadratic functions of v and its gradients. Thus when 
Vp is combined with the v ■ Vv term into dt"v, the third term of (8) contributes a triple 
correlation T^, similar to the triple correlation T^. in that it is composed of 2 powers of v 
and one power of s. We combine the two triple correlations and write T = T^ + T^. Deflning 
Q = —V • Vc, the coupled equations to be solved are then 

dtC^Q (10) 

^"""^ _ 1_ _ 1_ _ _ 

dtQ = gV'V^C + T = -v^V'C - Q/t,, (11) 

where in (11) we have replaced T by a damping term of the form —Q/ra, where = fTed 
and / is a dimensionless constant expected to be ~ 1. This replacement is our closure. 
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Note that since Q — — v • Vc = —V • (vc), for incompressible flows, Gauss' theorem then 
shows that Q represents a negative flux of s through the volume of the averaging. 

This "minimal r-approximation" for T and the value of / can be tested based on the 
implications we now derive. Combining (10) and (11), we have 

dJC + dtC/ra - ^^V^C = 0. (12) 

If we assume solutions of the form C = CoExp[nt + ik ■ x] then (12) becomes 
+ n/Td + l/rl = 0, where 1/tl = k^v"^ /\/3 ~ kiV2- The solution is 

" = ^(-i±\/i^)- 

In the limit that tl < Td/2, which corresponds to the case when triple correlations are 
assumed small, the solution becomes oscillatory. This would be unphysical if concentrations 
of passive scalars do not oscillate. To order of magnitude, this regime can also be written 
ki > ^2/2, which includes a regime in which the turbulent scale can be larger than the 
mean field scale, threatening the scale separation of the mean field formalism. 

We now consider the case tl > 2Td, or equivalently /ci < ^2/2. If /ci = A;2/2, 
then n ~ — l/(2rrf) as expected when ki and k2 are not widely separated. The 
more physical limit is fci << k2 so that » r^. In this case, we Taylor 

expand in Tj/r| to obtain n = ^ ^—1 ± ^1 — 7^^^- The "+" solution gives 

n ~ —Td/rl ~ ~{f /'^ed){k'l/k1) = {P / Tdijikl / k^) which gives the expected time scale of 
diffusion for a quantity of scale l/Zci, subject to turbulent motions of scale 1/A;2 when / ~ 1. 
The "-" solution, which is new, gives n ~ — l/r^. The total solution can be written 

C = Co,+Exp[-f(/VTd)(A;i/A;2)'] + Cc-Expf-i/rJ. (13) 

For / ~ 1, the first term on the right depends on the longer decay constant and will 
dominate for all finite 6*0,+ at late times. The second term on the right has a much 
shorter decay constant since ki « k2 and requires further analysis. Toward this end, we 
solve for Co,_ and Co,+ for several choices of initial conditions for dtC\t=o — <5(0) and 
Co^C{t^'o)^Co,+ + Co,->0. 

First consider the case Cq > and dtC\t=o = Q(0) = 0. From (13) we then obtain 
Co,+ = -{k2/kifCo,- and thus Co,+ = Co/{l - kj/k^). Then 

C - ^i_%ik2f M-tUVra){k.lk2f] - (rrf^Exp[-(t/r,)]. (14) 
The first term on the right dominates for all time, and the solution is thus similar to the 
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familiar diffusion solution of (4). Prom (10), the time derivative of this solution gives 

9.C = Q = -gig|^Exp[-t(/7r.)(fc,/fc,)T + - «ZM)^Exp[-(t/rJl. (15) 

At early times, the flux Q will begin to deviate from zero. In particular, the second term on 
the right decays more rapidly than the first so Q evolves to be negative. This corresponds 
to a flux of s from outside to inside the volume over which the averaging is performed. 

Let us now consider the case in which Q is initially finite, meaning that there is an 
initial net fiux of c in or out of the averaging volume. More specifically, we consider Co = 
and dtC\t=o = Q{0) 7^ 0. We then obtain Co,_ = -Co,+ = -TdQ{Q)/{l - kl/kl). Then 

C - j0§^^M-t{f/ra){h/hf] - (l?H|yExp[-t/r.], (16) 

again we see that C diffuses in the usual way; the first term on the right dominates for all 
time just as in the previous case. If Q(0) > ((^(0) < 0), then C becomes increasingly less 
positive (negative) with time, asymptoting to zero. Taking the time derivative gives 

^ {i^^WW /^d){ki/k2) ] + ^i_j,2j}^2^ ^M-t/rd]. (17) 

For t < -2Td'Ln[{f{kJk2)]/[l - P{kl/kl)] and k^/k^ « 1, (17) becomes 
Q ~ Q(0)Exp[-t/Td]. For Q(0) > (Q(0) < 0) the initial net fiux of s is outward 
(inward) and the solution thus implies that this outward (inward) fiux would decay on a 
time scale ~ before the first term on the right of (17) takes over. Then Q changes sign, 
becomes negative (positive) and eventually diffuses on the same time scale that C diffuses. 

For / ~ 1, it is noteworthy that in the above two sets of solutions for different initial 
conditions, C essentially diffuses as it would if it it were governed by (4), the "textbook" 
equation for mean scalar diffusion. But the evolution of Q in our approach reveals more 
subtle properties of the transport process. For an initially finite Q and Co = 0, the fact that 
Q incurs an initial decay in the presence of isotropic turbulence on a time scale ~ T^a can 
be understood physically: Without an initial mean field such as C, no net flux of s can be 
expected to survive longer than an eddy turnover time for isotropic nonhelical turbulence 
since the assumption of statistical isotropy would be otherwise violated. But as C becomes 
positive (negative), Q evolves to represent an inward (outward) flux of s associated with 
the diffusion of C and then evolves on the same time scale as C. 

Further insight into the closure T — Qjr^ comes from thinking qualitatively about 
entropy. A state with C = but large magnitude of Q can be thought of as one with 
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lower entropy than one in which Q has a small magnitude, since a large Q represents a 
macroscopically ordered anisotropic state. Increasing the entropy with time would mean 
decaying Q. This motivates a closure in which T is replaced by a term representing decay 
of Q rather than growth. The precise relation between T and Q can be tested numerically. 

Our closure is akin to the "r approximation" used by [19] for the anisotropic part of 
turbulent quantities. We call our closure the "minimal r approximation" as we use it for 
the full turbulent quantities. This closure is simpler than the Eddy Damped Quasi Normal 
Markovian (EDQNM) [7, 20, 21] and the Markovian Random Couphng Model (MRCM) 
[22] closures. In these closures, dtT oc Q rather than T cc Q. In the MRCM however, like 
ours, the decorrelation of triple moments is taken to be time independent. For purposes of 
calculating the diffusion of a mean scalar, we do not worry about the effect on the energy 
spectral index resulting from this choice. Further numerical work should compare the 
relative effectiveness of these various approaches. 

Three other features of our approach are important to emphasize. First, we did not 
need to introduce a correlation time as correlations between quantities at different times 
never enter the theory. Second, the oscillations in the C found when triple correlations were 
assumed to vanish highlight an likely unphysical consequence of any mean scalar transport 
theory which ignores terms providing the turbulent cascade. The latter is provided by finite 
triple correlations, which then also suppress these oscillations. Third, with respect to the 
time evolution of Q, the present approach is more useful than that implied by (5). There, 
for Cq = 0, nothing happens without an additional source term, even if dtC\t=o = Q{0) 7^ 0. 
Thus, unlike the system of (10) and (11), (4) does not model the evolution of a system in 
which a concentration of material is fed into a turbulent flow initially devoid of material. 

Our approach herein has now been shown to be consistent with fully 3-D numerical 
simulations of turbulent diffusion of a mean scalar (24). But the approach is also applicable 
to any equation for a mean quantity that depends on the correlation of two fluctuating 
quantities. A generic summary of the procedure is to replace the flrst order equation of 
the quantity whose time derivative is of interest by a second order equation by taking the 
derivative. The derivative of the double correlation is then to be re-expressed using the 
known dynamical equations for the time derivatives of each of the factors in the correlation, 
and then one takes the triple correlations that arise to be linearly proportional to double 
correlations. The coupled system of equations must then be solved. The approach has 
also been successful [16] in mean fleld electrodynamics, where the time evolution for the 
mean magnetic helicity depends on the curl of the turbulent electromotive force £ = v x b, 
where b is the fluctuating magnetic fleld in Alfven units. The procedure involves taking 
the time-derivative of £, using the equation of motion for dtV and the induction equation 
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for dth, and replacing the resulting triple correlations by a damping term oc £. The result, 
using Td ~ Ted if ~ 1) is consistent [16] with MHD simulations [23]. 

Thanks to A. Brandenburg for comments. EB acknowledges DOE grant DE-FG02- 
00ER54600. 
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